Spatial entanglement and massive neutrino oscillations produced by orbital electron 

capture decay 
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The two-particle wave function of neutrino and recoil nucleus is found as a solution of an initial 
value problem in the far zone for a time longer than the electron capture decay lifetime of a hy- 
drogenlike ion. The neutrino-recoil entanglement arising in such a process is a consequence of the 
momentum conservation and is closely related to the wave packet structure of the state. Because 
of neutrino mixing, the joint wave packet involves the coherent superposition of the neutrino mass 
eigenstate packets. This is the new physical realization of the Einstein-Podolsky-Rosen thought 
experiment, which has no analogue in quantum optics and quantum informatics. A class of possible 
experiments for the registration of a neutrino and a recoil nucleus is proposed. It is shown that due 
to spatial correlations neutrino oscillations can be observed in the coincidence experiment with the 
recoil. 

PACS numbers: 14.60.Pq, 23.40.-s, 03.65.Ud 



I. INTRODUCTION 

Massive neutrino oscillations are a consequence of the 
presence of flavor neutrino mixing and a clear evidence 
of physics beyond the standard model. They are sub- 
ject of an intense experimental and theoretical research 
beginning with the pioneering paper by Pontecorvo [ij. 
The theoretical approach is mostly phenomenological and 
uses the plane wave 0-3 or the wave packet formalism 
d-Q to describe the evolution of the massive neutrino 
states, and the field-theoretical approach flo'-^T^ which 
takes into account the processes of production, propaga- 
tion and detection of neutrino. In spite of the fact that 
almost all these treatments provide the canonical formula 
for the probability of oscillations, some basic issues of 
the theory of neutrino oscillations are still being debated 
[isj . In particular, there is no consensus on whether the 
three neutrino mass eigenstates have equal energies or 
equal momenta. Furthermore, the wave packet formal- 
ism leaves unanswered the question of how the properties 
of a neutrino wave packet are determined by the process 
of weak decay. 

If we want to resolve these problems, we need to do 
away with ad-hoc assumptions and invoke basic princi- 
ples of quantum mechanics to describe the weak decay of 
an unstable object and neutrino production • Closely 
related problems have been studied in quantum optics. 
Fedorov et al. [l^, [11] considered the decay of composite 
objects into two fragments that are free to move away 
from the breakup point and are constrained only by mo- 
mentum and energy conservation. The authors obtained 
position dependent two particle wave function as the so- 
lution of an initial-value problem. This quantum state is 
entangled and closely related to the Einstein-Podolsky- 
Rosen one 



'Electronic address: ' paviQmbslab.kiae.ru| 



Our motivation for this work is to apply the formal- 
ism of Refs. [l^ [iBl to the description of orbital electron 
capture (EC) decay. In this scenario, the spontaneous 
emission of a flavor neutrino by an atom is described in 
three dimensions with initial wave function of a decay- 
ing atom taken in the form of a finite-size wave packet. 
The recoil-neutrino wave function is found analytically in 
the coordinate representation as the solution of the time- 
dependent Schrodinger equation in a far zone. Its wave 
packet structure evolves with time. The function does 
not factorize in the neutrino, r^, and recoil, r„, spatial 
coordinates, which is a direct indication that the quan- 
tum state of the system is entangled. However, due to 
neutrino mixing the entanglement involves the coherent 
superposition of mass eigenstates with equal energies and 
different momenta that manifests itself in the oscillating 
correlations of the two particles. 

Entanglement means that knowledge of one of the par- 
ticles reveals information about the other. Space correla- 
tions of neutrino and recoil nucleus are the experimental 
implication of this feature. It should be emphasized that 
the spatial entanglement which we are considering, has 
nothing to do with kinematic entanglement (see, for ex- 
ample, Refs. 0, Hj) which does not involve an exact so- 
lution of the Schrodinger equation. Spatial entanglement 
between an atom and a photon in spontaneous emission 
has been observed by Kurtsiefer et al. [l^ . 

Because of strong spatial correlations, we are able to 
observe with two detectors neutrino oscillations and a 
recoil simultaneously in the coincidence measurement. 
This experiment, first considered by Dolgov et al. [20j . 
could be of interest for the next generation of neutrino 
detectors. By contrast, registration of one particle in- 
dependent of the other disentangles the recoil-neutrino 
pair. For example, to observe neutrino oscillations the 
recoil position is ignored whereas a flavor neutrino is de- 
tected. This scenario is used currently to observe neu- 
trino oscillations in reactor and accelerator experiments. 

The paper is organized as follows. In Sec. II we ap- 
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ply the Weisskopf-Wigner theory of spontaneous emission 
|21l | to treat the EC-decay process of the ion confined to a 
small volume in a trap. We specifically study a hydrogen- 
like (H-like) ion of the intermediate mass region with a 
K-electron in the lowest hyperfine state having the total 
angular momentum = / ± 1/2, where / is the nuclear 
spin. Such ions have been used in the much-talked-of GSI 
experiment [22| and its second run in 2010 [23]. However, 
the question of what happens to the neutrino and the re- 
coil nucleus on time scales comparable to the lifetime of 
the parent ion is left unanswered. The main result of this 
section is the recoil-neutrino wave function that describes 
the spatiotemporal behavior of these particles. In Sec. Ill 
we use the far zone approximations to express this func- 
tion as a product of the relative and the center of mass 
(CM) wave functions. It is shown that the coherent su- 
perposition of the neutrino mass eigenstates has a fixed 
kinetic energy equal to the EC-decay energy. The struc- 
ture and time evolution of the relative motion (RM) and 
the CM wave packets are investigated in Sec. IV. The far 
zone formalism allows to find the spreading of these wave 
packets. In Sec. V we outline the experiments appropri- 
ate for the detection of a flavor neutrino and a recoil 
nucleus. Our findings are summarized in Sec. VI. The 
Appendix contains the details of calculation of an inte- 
gral from Sec. III. 



II. ORBITAL ELECTRON CAPTURE DECAY 
OF MOVING HYDROGEN-LIKE ION 

The decay which we study is the capture of a sin- 
gle electron from the K-shell into a completely ion- 
ized daughter nucleus and a electron neutrino, as fi- 
nal state. The system is described by the Hamiltonian 
H = Hq + Vsf + Hyj , where the unperturbed Hamiltonian 
has the form 



P 

2M 



Hn{pi,...,pA)+^exa 



tax 



ei(fc)c+Cji; 



(1) 

Here p is the nucleus momentum operator; is the 
nuclear hamiltonian depending on nucleon coordinates, 
Pk = Tfe — r„, with respect to r„ which to a good approx- 
imation is the nucleus CM; A is the mass number; 
and a\ are the creation and the annihilation operators of 
a bound electron in the state with the energy e\ and the 
quantum numbers X — n, j, m; Vsf is a hyperfine interac- 
tion; and and c^k are the creation and the annihila- 
tion operators of the three massive neutrinos {i = 1,2,3) 
with momentum k, the energy = y^k'^ + mf and the 
mass rrii. The mass M of parent and daughter particles 
is supposed to be equal with an accuracy of the small 
parameter Qec/M, where Qec is the decay energy. 
The weak-interaction Hamiltonian density is given by 



Gp 

7f 



Vud [j(r)-J(r)+H.c.: 



(2) 



where the lepton current, 

j(r) ^ ii?e(r)7(l +75)e(r). 



(3) 



involves the flavor neutrino and electron field operators. 
The former is [Hi 



.(r) = ^ UeiVe{v) = ^ C/eiCik.1tL(n,)e'''*' 



(4) 



where U is the Pontecorvo-Maki-Nakagava-Sakata mix- 
ing matrix. The Dirac spinor for left-handed neutrino, 

1 r-M...rn^.r-^"^^'^^-(^/2)V 



"^^") = 7!l"-J'""^"^"V e-/^cos(^/2) ; 

(5) 

depends on the unit vector n = k/fc. We will use non- 
relativistic limit for the electron operator 



(6) 



For the K-electron, we have Wj^i^ = fis{r)rii where 
/ is radial and 77 is spin wave functions. The nu- 
cleon current for pure Gamow- Teller transitions in a non- 
relativistic approximation has the form 

A A 

J{v)=gAj2 '^(^)^- (s)<5(r-r,) =.g^^ J(s)<5(r-r,), (7) 

s=l s=l 

where a and r are the spin and the isospin Pauli matrices. 

To find the wave packet structure of an entangled 
state, one needs to use the coordinate representation 
for the wave functions of neutrino and recoil nucleus. 
This implies that the Hamiltonian i?^, has to commute 
with the total momentum of a bipartite state in the r- 
representation. By using Eqs. ^ and (O we get af- 
ter integration over nuclear volume the weak interaction 
Hamiltonian fulfilling the momentum conservation law 

Hn. = z^KdJ] ;7*c+^UL(n)7(l + 75) 



: J(s)e-*''*''= aAV'A(p.)e"''('-"-'-") + H, 



c. 



(8) 



Now let us use the time-dependent perturbation the- 
ory to determine the temporal evolution of a decaying 
state for the Gamow- Teller transition FIMp I'M' 
with /' = /± 1. The parent ion is prepared in a polarized 
state with the angular momentum projection Mp along 
the quantized axis Z of the laboratory frame S. The so- 
lution of the Schrodinger equation with the Hamiltonian 
H is sought by using the following ansatz for the wave 
function 



= Y -^(P' ^) \aFIMp; p, t) i 



^ S,M'(p,k„0|W'M';p-k„k„t)e-'^^'*, (9) 
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where 



m.M 

\hl'M'- p - k„ k,) = ^M'M' (p) C+ . |0) e^(P-k.)r„+«k.r„ 



(10) 



are the eigenvectors of i/o + Vgf with the eigenvalues 

(P - 



2M 



(11) 



In Eq. ([TU)) ft{p) — il{pi, Pa) and E are eigenfunc- 
tions and eigenvalues of the initial (al) and final {bl') 
states of nucleus, Eq is the energy of the ion ground state, 
|0) is the lepton vacuum, and C'^2m /Af Clebsch- 
Gordan coefficients. In the subsequent text we will as- 
sume that Ef,!' = 0. 

The differential equations for the coefficients A and B 
are 



E 

i,ki,M' 



iB,M' (P, k„ t) = ^(p, t)WlM- (n»)C/* e 



where the matrix element of the Hamiltonian H,, 



,(12) 



F Vud 



[f is) gAMiaI,bI')C^A^i^^), (13) 



includes the radial wave function of the bound electron 
averaged over nuclear volume, (/is), and the nuclear re- 
duced matrix element Ai. The spinor 

6'M'(n)^C,^:^f,X'-i(^'^'0)' (14) 

where D is the Wigner function and /i — Mp^M' = ±i, 
is normalized by the condition X^a/^ m' \^I'M' \^ ^F+1/2. 

Suppose that a parent ion after production is stored in 
a trap and its motion is localized. Then we switch off the 
field of the trap and free spreading of the ion CM wave 
packet begins. If the production and stored time is far 
less than the life-time of the parent ion, we can measure 
time from the beginning of free spreading. The initial 
state is described by the first term of Eq. ^ ior t = Q. 
To be specific, suppose that the CM part of this function 
has in the momentum representation the Gaussian form 



A{v,t^Q) = {2^d)^/^ exp (-^rfV 



(15) 



The Eqs. ([T2|) are easily solved in the Weisskopf- Wigner 
approximation with the initial conditions ^(p,i = 0) = 
^o(p) and Sij\f'(p, k, i = 0) = 0. The solutions for t >• 
^IQec are 



yl(p,t) = A(p)exp(-R), 



Biw (p, k, t) = 
where 



Ao{p)W!,,,,{n,)U*^ 

£bi - £a+i^ 



(16) 



r = vr^ \WpM'{ni)f \Ue^f 6{£u - £a) 



i,M' Mi 



= ^^^^'^^MdAMial^bnQEcf (17) 

is one half of the rate of the Gamow- Teller transition 
FIMp I'M' , I' — I ±1. The final expression we have 
derived for this value ignores a small dependence of F 
on the parent ion velocity. Furthermore, the two small 
dimensionless parameters 



a 



Qec 
M '■ 



6^ 



Qec 



(18) 



allow to reduce F to a conventional value for the ion at 
rest. The last parameter is compatible with the limit of 
ultra relativistic neutrinos (25l |. 

First of all, we observe that the decay probability P{t) 
of the parent ion is not dependent on the shape of the 
initial CM wave packet ^o(p) because of 



P{t) = exp {^2Tt)J2^oiP) = exp (-2Ft). 



(19) 



However the probability of emission of the electron neu- 
trino with wave vector k and energy e is little affected by 
Doppler shift 



^e(k): 



Af 



\WrM'in)\ 



aUp) 



EC- 



2A/ A/ ) ' 



(20) 

With BiM' taken from Eq. p^ . the recoil-neutrino 
wave function at times t > 1/T has the form 



*/'A/'(r„,ry,t) = ^ 



^o(p)W^/V(n^ 



M 



j,p,ki ^i — QeC- 

X exp{i(p-ki)r„+ ikiVu~i£bit} Ul^ct^, |0) 



(21) 



It carries information on the neutrino production process 
and involves the coherent superposition of the electron 
neutrino mass eigenstate components. 



III. RECOIL-NEUTRINO WAVE FUNCTION 

Now, by using the far zone approximation and small 
parameters (|18p . we write the function (j2ip in the form 
suitable for the analysis of the wave packet structure and 
entanglement of the recoil-neutrino state. To begin with, 
let us rewrite Eq. (PT|) as 

"Hi'M' =^A(p) exp (*Prn-«|^] E ^) ^e»cijO) , 

p ^ i 

(22) 
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where the function 

> WIm' (n» ) exp ik,; r - i (e, + ij. - 2^ 



Xi 



' 2M " 



M 



(23) 



depends on the relative coordinate r = — r„. In the 
laboratory frame S*, the vector r is specified by polar 
angles 9 and ip. The direction of ki is given by angles "di 
and (pi. 

First we transform the summation over kj in Eq. 
to integration. To perform integration over dTl-^^ , we ro- 
tate the laboratory system S with the origin O in the 
decay point by the angles 0, Lp to align axis Z along the 
vector r. In the new frame S", the vectors k^ and p are 
determined by angles 1?^, 0^ and ??p, 0p, respectively. The 
polar angles of these vectors in the systems S and S' are 
connected with each other by well-known formulas [26| . 
Owing to kiT ^ QEct > Qec/^ ^ 1, the integrant in- 
volves rapidly oscillating function exp (i/cir cos t^^) and a 
slowly varying preexponential function T of the angle 
Integration by parts over this angle yields 



F{-d[) exp {ihrcos^'i) sin<d< 



—iki 



J-(0)e''=''-] +0(l/(fc,r)2). (24) 



Here the two terms in brackets represent incoming and 
outgoing spherical waves. The incoming wave in the 
far zone gives an exponentially small contribution for 
t > 1/r which can be neglected in comparison with the 
outgoing one. The latter corresponds to i?- — 0, and it 
immediately follows that di = 9 and (pi = p 'm the inte- 
grant of Thus, we find that ki||r and = n = r/r 
in the far zone. After integration over d(j)'i, one obtains 



exp 



ikiV^i [ ei + 



2M 



M 



-Qec- 



' 2M M 



iT 



(25) 

where pki = pki cos -dp. 

To calculate the integral, we change variable ki by 
and expand the former around Qec 



ki{ei) = fcjoH — i^i -Qec) - 



{er-QEcf, (26) 



because main contribution to the pole in (j25p comes from 
— Qec — ir. We have retained a quadratic term to 
describe the spreading of the RM wave packet. Here 

ho = V Qec ~ "^f ^'^'^ ^» = ho /Qec is the group ve- 
locity of neutrino. Then the integral over can be eval- 
uated by the residue method. The details of calculations 
are given in the Appendix. Upon using Eqs. (|13|) and 
(|17p , we find in the lowest order in the parameters a and 
6,. 



Xi 



-ipi{n, r,t) exp(iapr). 



(27) 



where the RM wave function of ith massive neutrino is 
~1 



V'z(n,r,i): 



1-Erf 



if. Vit — r 
2 V^"^' AZd" 



X ^i,M' (n) exp [iihar - Qsct) - [vit - r)/A]- (28) 

Here Erf is the error function, Di — Vi/T is the initial 
width of the RM wave packet, and 



1 



t 

M 



_6fr_ 

rriiVi 



1/2 



(29) 



is a dimensionless parameter describing its spreading. 

Now, integration over dp can be readily performed af- 
ter substituting expressions ([T5|) and (P7|) into Eq. 



E-^o(p) 



exp 



zp r„ +av]~i 



r3/4 



3/2 



exp 



(r,i 



2d(d+^ 



(30) 



In view of the definition of the CM vector for a relativistic 
system in terms of the particle energies [l^ , 



R 



(31) 



we find an analytic expression for the wave function (j22p 



*/'M'(R,r,i) = 



1 



3/2 



exp 



R2 



^V.(n,r,i)t/*c+o|0) 



(32) 



The two-particle wave function (|32p describes the evo- 
lution of the recoil-neutrino state after decay at time 
t > 1/r. This expression is accurate to within the small 
parameters a and Si (jlSp . The function carries informa- 
tion on the decay process and has the form of a product 
of the CM and RM parts. Such a factorization of the to- 
tal wave function is a general feature of the two-particle 
decay of noninteracting fragments constrained only by 
momentum and energy conservation [T6| . The RM wave 
function involves the coherent superposition of the mass 
eigenstate components of the electron neutrino state. 

The time-dependent phase factor of the recoil-neutrino 
wave function ([5^ [see also Eq. (^5]) ] involves the kinetic 
energy of the pair. As shown in the Appendix, it is equal 
with great accuracy to the decay energy Qec- Thus, 
all three massive neutrinos have the same energy. This 
resolves the long standing problem " same energy or same 
momentum." This issue, like other paradoxes of neutrino 
oscillations emerges from the theory that considers 
the time evolution of a neutrino independently from a 
recoil. The exact solution (15^ shows that a neutrino 
and a recoil do not evolve separately due to their spatial 
correlation. 
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IV. THE WAVE PACKET STRUCTURE OF 
THE RECOIL-NEUTRINO STATE AND SPATIAL 
ENTANGLEMENT 

The spatiotemporal behavior of the joint quantum 
state of the recoil and the electron neutrino following 
the EC-decay is in agreement with the results obtained 
in the theoretical studies of decaying bipartite systems 
[isl [l6l , \2l\ . The distinctive feature of our system is the 
coherent superposition of neutrino mass eigenstates. In 
this section we give the more detailed analysis of the 
function p2p providing insight into the nature of entan- 
glement and neutrino oscillations. We begin with the 
wave packet structure of this function, since entangle- 
ment and neutrino oscillations depend on spatial local- 
ization of particles involved in the decay. 

The CM part of ^P/'M' (R-, r, t) has the form of a spread- 
ing wave packet 



*CM(R,i) 



1 



7r3/2D3(t) 
with the time-dependent width 



exp 



R2 



(33) 



(- 



1/2 



/2 



d, t < tcM 
/d, t > tcM 



(34) 



where tcM — Md? is its spreading time and Im — \/t/M 
is the quantum diffusion length. If the initial size d of the 
atomic wave packet is approximately 10~^ cm, we have 
tcM ~ 10~^ s for nucleus with the mass number A ~ 100. 
Hence, the width is of pure dispersion origin and 
grows linearly with time. The velocity of spreading is 
equal to l/Md. 

The time-dependent width of the RM wave packet 
|V'i(r, i)P is due to dispersive broadening of both the re- 
coil and the massive neutrino. To estimate these effects, 
let us put r « Vit in Eq. ([29]). We get 



A,: 



r3,2 Nl/2 hl_ 



(35) 



where ti — MD^ is the spreading time caused only by 
a recoil particle. The contribution of a massive neutrino 
is negligible because of the relativistic suppression of the 
wave packet spreading. For F ~ ls~^ and A ~ 100, the 
spreading time is U ~ 10^^ yr. This result is apparent 
for ultrarelativistic neutrinos, and we can put Ai = in 
Eq. (I^S)) . In this limit, one finds the RM function 



V'/'M' (n, r, t) - 



1 



/2ti 



e/'M'(n)^i?.(r,i)f/*c+^jO), (36) 



where the radial function of the ith massive neutrino is 



. X v2F 
Ri(r,t) = exp 



iihor-QEct)- 



Vjt-r 
D, 



(37) 



Here Q is the unit step function. The function Ri (r, t) is 
Lorentz covariant. The RM function p6p is normalized 
for ultra relativistic neutrinos by the condition 



[drY,\ i^i'M'{n,r,t) ^ 

M' i 



1. 



(38) 



The wave packet | ipi'M'{^,t) P is a superposition of 
three exponential wave packets of massive neutrinos with 
different sharp edges r — Vit and widths Di. The later 
depends only on the dynamics of the decay process. The 
difference between the group velocities Vi of these pack- 
ets results in their separation. However, the separation 
is negligible compared with Di for times t ^ lO^'^ yr. 
Therefore the wave packets corresponding to different 
mass eigenstates are in fact spatially inseparable. Hence, 
we can use an ultra relativistic approximation for Ri with 
Di — D = 1/V and Vi = 1. In this approximation the 
function ([57]) takes the form 



Ri{r, t)^Ro{r, t) cxp[i(fcior-(3£;ci)]: (39) 



where 



i?o(r,i)- Jy|exp(^-y^je(i-r). (40) 

We now wish to consider the two-particle wave func- 
tion, taken as the product of two functions (|33|) and p6p , 
in the observable coordinates of the recoil nucleus and the 
neutrino. This normalized function has the form 



I'M' 



CM 



((l-a)r„4ar^,i)-i/'/'Af'(n, |r^-r„|,t). (41) 



The function does not factorize in these variables - a 
direct indication of the spatial entanglement of two par- 
ticles. Each of the three massive neutrinos becomes en- 
tangled with the recoil nucleus because a neutrino is not 
emitted in a momentum eigenstate. 

The joint recoil-neutrino wave packet X^m'I^^™'!'^ 
has an axially symmetrical shape with respect to the axis 
passing through the CM in the direction of the vector 
n. The packet increases with time in a transverse di- 
rection with velocity l/Md and in a longitudinal (along 
the axis) one with velocity Un + (1 — w„) = 1, where 
Vn = Vi^i/^n ~ ctVi Is the velocity of a recoil. It is easy 
to see that the function ^cm takes its maximum value 
on the symmetry axis, along which the state is highly 
entangled. The probability density is proportional to a 
product of the Gaussian and exponential functions 



El* 



I'M' I ^ Fp {9} exp ■■' 



M' 



xexp[- ^ '"" je(t-r.-r„), (42) 

where the angular modulation of the joint packet is de- 
termined by the function (angular distribution function) 



^ ' 47r I /' + 1 



for/'=/±l. (43) 
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Modulation is due to the polarization of a parent ion. 
For a non-polarized ion we have Fji — l/47r. 

The Gaussian packet in Eq. ([^^ for fixed has the 
form 



exp 



1 - a 



whereas for fixed r„ it is equal 



exp 



D2 



(44) 



(45) 



The probability distributions (|48)) and p9)) reveal no en- 
tanglement because all information about the position of 
one of the particles is lost completely. 

Coincidence measurement. On the axis with highest 
entanglement, the probability density to detect recoil nu- 
cleus at the point r„ together with the neutrino of flavor 
P at the point is 



dP, 



e/3 



drndr„ 



= ^^/'(^)*cm((1 



a)r„ -I- ar^,t)rlrl 



^ Up.U:,Ue,U;^R* (r. + r„)i?, (r, + r„). (50) 



Relative location of the peaks of these curves is deter- 
mined by the condition 



(1 - a)r„ - ari, = 0, 



which corresponds to the maximum of the wave packet 
(P^ . The dependent exponential and Gaussian curves 
overlap each other, if r„ does not exceed its maximum 
value Vnt. 



V. EXPERIMENTAL IMPLEMENTATIONS 

Recent experiments concerned with atom-photon en- 
tanglement are dealing with generation and verification 
of an entangled pair. The special structure of the recoil- 
neutrino wave packet offers new kinds of experiments. 
There are three types of possible experiments to detect 
a flavor neutrino and a recoil: 

(i) Coincidence measurements, in which both the re- 
coil and neutrino are detected. The flavor-changing 
process e — >■ /3 is determined by the probability den- 
sity 



(0| ^Cji,^oUpj'^i'M'{rn,r^,t) 



3 

M' 3 



(47) 



(ii) Noncoincidence measurements, when only neutrino 
is detected regardless of the recoil position. The 
flavor-changing process e — /3 is determined by the 
probability density 



/ dr„ ^|(0| ^Cji^.„Ui3j'^rM'{rn,ri,,t) 

M' j 



(48) 



In ultrarelativistic approximation the second line is the 
probability of the e — >■ /? transition 



(46) E \Ue^mp^\'+2Y,\up,u:,u,,u*p^\ cosi 



i>j 



2Qec 



r+7 



(51) 



where Ainf, = m: 

''J 



ni] and 7 = a.Tg{UpiU*^UejUL 



We choose the direction of the axis corresponding to 
the maximum value of Fp(0), which will be denoted by 
Fj'. The neutrino detector N, which fixes the flavor of 
neutrino states by some charged current process, and the 
recoil one R to register the nucleus are connected with a 
coincidence circuit and located on this axis on both sides 
of a parent ion confining volume (that is approximately 
the CM position) at distances from it, respectively, L^, 
and Ln- If and L„ satisfy Eq. (j46| . the amplitude of 
oscillations will be equal 



TT^I'^DD\{t) 



exp[-2r(i -Ly-Lr. 



(52) 



The amplitude conforms to highest spatial entanglement 
and can serve as its verification. The detection time t 
depends on the lifetime of a parent ion. For short life- 
times, the length of the joint wave packet may be equal 
to distance between the detectors, that is i = L^+L„. 
In this case, the amplitude (j52p will be maximum. This 
condition is impossible for parent ions in a Penning trap 
because their lifetime must exceed the time it takes 
to prepare them in the trap. The latter is of the or- 
der of seconds. In such a case, the measurement time 
t> XjV should be taken so as to minimize the exponent 
in Eq. ([5^ . 

It follows from Eq. ([5T|) that the expression for the 
probability (|50l) contains the well-known oscillation phase 



(iii) Noncoincidence measurements, when only the re- 
coil nucleus is detected regardless of the neutrino 
position. The probability density to detect recoil 
nucleus at the point r„ is 



'' M' j 



(49) 



2Qec 



For two neutrino mixing, i.e. in the case of e 
/3 = e, /i) oscillations, the mixing matrix is 



U = 



cos rj sm rj 
— sin 77 cos 77 



(53) 



(i.e. 



(54) 
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where rj is the mixing angle. Then, one finds from 
Eq. ([5T|) the probabilities 



Wee - l-sin2(2r/)sin20i2, We 



sin^(277) sin^(/)i2, 



(55) 

where (/)i2 = 7r(ii/ + i„)/Lo and Lq = ^t^Q ec I ^"^"^ 
(Am^ = m| — > 0) is the oscillation length. The 
final expression for the coincidence probability density is 



e/3 



2FrLl 



X exp[-2r(i - i,-L„)]We/jW„, (56) 

where W„ = 1 is the probability of a recoil registra- 
tion. For the detectors connected to coincidence circuit 
with one output (the "and circuit"), a signal disappear- 
ance at certain positions of the detectors N and R is 
caused by oscillating probability We^j, in the muon neu- 
trino channel. The recoil nucleus does not oscillate. To 
observe the spatial correlation and neutrino oscillations, 
one needs the correlation experiment, in which an exper- 
imental event involves the registration of a neutrino si- 
multaneously with a recoil nucleus in two detectors. One 
can fix the distance from source to one of the detectors 
and change the distance to another one. Here, too, both 
neutrino oscillations and a recoil may be observed, but 
with a far lower amplitude due to the Gaussian factor of 
Eq. gSl). 

Let us examine the conditions for the observation of 
oscillation patterns. The first condition is trivial: the 
source-detector distance should be of the order of, or 
greater than, the oscillation lengths Lose- The second 
condition is the coherence of different mass eigenstates, 
necessary for the neutrino and recoil oscillations to be 
observed. Coherence is preserved over distances not ex- 
ceeding the coherence length Lcoh- The latter is defined 
as the distance at which the phase difference due to en- 
ergy spreading obeys the equation 

<^>^J (Qec) - c^^J {Qec + T) = 27r. (57) 

We find Lcoh = LoscQec/^ ^ lO^^iosc- Such a large 
correlation length arises because the emitter is in a pure 
quantum state that is described by a state vector ([9]). 

Neutrino detection. This is a well-established experi- 
ment for studying neutrino oscillation. Integration over 
dr„ in Eq. (I48p can be done approximately since the 
width of the Gaussian packet p4|) is much smaller than 
that of the exponential, if Ft <C 10^". Indeed, we have 



t^ 



Ft 

{l~a)D Z? V " ' JPd? 2Md 
The integral is calculated as follows 



10"i"rt. (58) 



*cm((1 - a)rn + ar^)R*{\r^ ~ r„|)i?j(|r^ - r„|)dr„ 



R* 



Let the position of a neutrino detector be given by the 
distance from the emitter, then the probability den- 
sity to detect the neutrino of flavor f3 at this distance 
is 



dPeB 2Fi 



'-el^v[^2T{t-L,)]\J2\Ue^\''\Up 



-2Y}Up,U*,,Ue,U*p^\cos(2^^ 



+7 



(60) 



Recoil nucleus detection. Integral (I49|) is calculated in 
perfect analogy to P5|) if Ft <C 10^ As a result, we obtain 
the following expression for the probability density to 
detect recoil at a distance L„ — Vnt from the parent ion 



dP, 



FrT 



drndVln 



(61) 



The joint wave function of the recoil-neutrino pair al- 
lows to calculate probability densities of the considered 
processes. For full description of oscillation experiments, 
it is necessary to take into account a detection process. 
However, this is different problem. 



VI. CONCLUSION 

In summary, we have found an accurate analytical so- 
lution for the joint quantum state of an electron neutrino 
and a recoil nucleus following the electron capture decay 
of a hydrogenlike ion. The evolution of the state provides 
an exactly calculable illustration of the famous Einstein- 
Podolsky-Rosen thought experiment. The new effect is 
the entanglement between the recoil and the coherent su- 
perposition of the three massive neutrinos. Such a super- 
position cannot be observed in the experiments involving 
discrete states of an atom-photon pair because the polar- 
ization of an emitted photon is inevitably entangled with 
the spin state of an atom. 

We have shown that each of the three massive neu- 
trinos is not emitted in a momentum eigenstate, and so 
the recoil becomes entangled with the superposition of 
three neutrinos. However, the total energy is conserved 
in the course of EC decay, and the decay energy converts 
to the kinetic energy of the reoil-neutrino pair. Thus, we 
have rigorously proved that the neutrino mass eigenstates 
composing the electron eigenstate produced in EC decay 
have the same energy. It should be emphasized that al- 
though a particular type of decay has been treated, the 
results of our calculations are applicable to other two- 
body weak decays due to the common features of the 
initial- value problem for nonintcracting particles. 

The most peculiar aspect of our solution is the two- 
particle system in a pure quantum state, where the par- 
ticles, a recoil and a massive neutrino, are mixed with 
each other. Such an entanglement opens up an oppor- 
tunity for correlative experiments and the investigation 
of environment-induced decoherence. This is the main 
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distinction of our approach from the field-theoretical one 
[Til] , where a neutrino enters as unobserved intermedi- 
ate state. We have suggested the correlative experiment 
which allows to observe neutrino oscillations and a re- 
coil simultaneously. However, a considerable progress in 
the detection methods of neutrinos and recoil nuclei is 
necessary to carry out such measurements with a single 
ion. 



Acknowledgments 

The author is grateful to A. L. Barabanov and Y. Litvi- 
nov for critical comments. This research was supported 
by Grant NS-7235-2010.2 from the Russian Ministry of 
Education and Science. 



Appendix 



Then, we first evaluate the integral 
^ ki{e) expliVie) + i{e - Qec)x] 
c 



de. 



(A.5) 



Vi{e){e~ep) 

where the contour C encloses the pole (|A.2|) . and the 
exponent V is 



P(e) = hor - Qec 



2M 



pn 



Mvi Mv, 



M 



{e-Q 



EC) 



(A.6) 



The integral is easily calculated, and Eq. (jA.ip takes the 
form 



I: 



2i7rfcjoe 



r jvi—t 



exp 



Vix{tp-QEc) 



dx, (A.7) 



We consider the integral 
"° h{e) exp|ifc,(e)r- 



I 



2M 



Pk.(0 

M 



fc2(e) pk.(£) 



-de. 



2M 



M 



(A.l) 



where the function ki{e) is given by power series ([26]). 
In the lowest orders of the small parameters a and for 
T ^ Qec the pole of the integrand is 



f^p = QEC-i^ 



( pkii 



V M 



2M J\ 



pn 



ho 



Mvi Mv, 



(A.2) 



We may extend the lower limit in (lA.ll) to — cx), since the 
contribution to the integral falls of sharply with increas- 
ing \£i\ owing to r <C Qec- The integral involves the 
exponential factor exp[—iqi{€ — QEc)^] with 



*=2 



6f 



Mv. 



1 / t 



M M 



5h 



2 \Mv^ rriiV- 



(A.3) 



which prevents the use of the residue method. To over- 
come the problem, we have to use the Fourier transfor- 
mation of this factor UM 



exp [—iqe 



exp 



4^ 



iex j dx. 



(A.4) 



where for F <^ Qec the expression for the exponent 
taken in the pole has the form 

^-ne,)-k,Ai--+-\-QEc{i+^-^At 



'da „ ^Vit 

apr 1 h2a2 — 

2 r 



(l-a^)t (1-a) 



+ a 



(pn)2 

M 



r 5a \ {pnf 

7rt + 01——:^ — t. 

V, 2 / M2w,, 



We have approximately 

V w kior + apr - Qsct + iT{t ~ r/vi). 



(A., 



(A.9) 



The last two terms in the expression (|A.8[) can be 
dropped because they are small compared to the term 
proportional to p^ in the exponent of Eq. . It is seen 
from Eq. (jA.Sp that the energy of a pair of ith neutrino 
and a recoil is the same for all pairs with great accuracy 
and is equal to the decay energy Qec- Small deviations 
from Qec are attributable to the approximate calcula- 
tion of the pole (jA.2p : the more precisely it is defined, 
the less are these deviations. 

Finally, the integral (|A.7p can be transformed into the 
error function 



inkio 



1 _ E^r-M-^zg^F-^ 
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